Nucleon-nucleon potential in finite nuclei 
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We consider the spin-isospin-independent central part of the residual nucleon-nucleon potential 
in finite spherical nuclei taking into account the deformation effects of the nucleons within the 
surrounding nuclear environment. It is shown that inside the nucleus the short-range repulsive 
contribution of the potential is increased and the intermediate attraction is decreased. We identify 
the growth of the radial component of the spin-isospin independent short-range part of the in- 
medium nucleon-nucleon interaction as the responsible agent that prevents the radial collapse of the 
nucleus. 
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I. INTRODUCTION 

Nearly 20 years ago, in their pioneering work m , Jack- 
son, Jackson and Pasquier investigated the skyrmion- 
skyrmion interaction by considering simple deformations 
of the chiral field. This idea was further developed in 
Refs. H, H, Q into various versions of Skyrme-like mod- 
els allowing for more complicated deformations. It was 
shown that the central potential of the original Skyrme 
model is reduced by an amount of about 10 20 % if de- 
formation effects of the skyrmion are taken into account. 
In the w-stabilized Skyrme model the deformation effect 
is even more significant, as the central repulsion at i? = 
Ifm acquires a reduction of about 40 %. 

Nowadays, in-medium properties of the nucleon are 
of high relevance. In particular, properties of a sin- 
gle nucleon embedded in various finite nuclei have been 
considered in a medium-modified version of the Skyrme 
model in IS' It ha-s been shown that bound nucleons 
acquire an intrinsic quadrupole moment due to the in- 
medium deformation effects. In summary, the general 
behaviour of the properties of the nucleons are in qual- 
itatively agreement with experimental observations, e.g. 
the swelling of the nucleon in the nuclear medium and the 
decrease of its mass. Similar results have also been ob- 
tained in the infinite nuclear matter approach of Ref . ■ 

Thus, it is only natural to investigate the nucleon- 
nucleon (NN) interactions in the nuclear medium because 
of their role in the formation of nuclear matter. There 
are two known paths to nuclear matter aspects of the 
Skyrme model. The first one is based on studies of the 



to studies of Skw'me model on the hypersphere 0, 0, 
El m mm 123. Especially in the first approach, the 
tensor part of the potential is found to be responsible for 
the crystalline structure of the ground state. In contrast 



to this, quantum hadrodynamics studies |24| assert the 
important role of the spin-isospin independent scalar part 
of the NN interactions in the formation of the nuclear 
matter. In any case, one can expect that the surrounding 
nuclear environment will leave its mark on the in- medium 
NN potential. 

The temperature and density dependence of nucleon- 
nucleon interactions, hadron properties and meson- 
nucleon coupling constants of the one-boson-exchange 
potential have been studied in the framework of thermo 
field dynamics (TFD) in Refs. HI El mm El. In par- 
ticular, it has been shown that the potential well of the 
NN interaction becomes shallow as the temperature or 
density increases 12^ . The changes in the meson-nucleon 
couplings and hadron properties also indicate the influ- 
ence of the medium on the one boson exchange poten- 
tial j28Ll29| |. Note, however, that the Skyrme model and 
its variants are the only type of models that use hadronic 
degrees of freedom and that allow for a simultaneous de- 
scription of both single-baryon and multi-baryon proper- 
ties. In fact, as discussed in Ref. |5i|, Skyrme-type models 
even allow for a simple inclusion of nuclear background 
terms. 

The alterations of the tensor part of the NN potential 
for infinite nuclear matter have already been studied in 
the framework of the in-medium modified Skyrme model 
in Ref. (Tj. In the present paper, we rather concentrate on 
the spin-isospin-independent part of the NN potential in 
finite nuclei where we allow the chiral fields of the two in- 
volved nucleons to deform in response to the surrounding 
nuclear environment. 

For simplicity, we will only consider the residual in- 
medium NN potential in the radial direction of the nu- 
cleus (i.e. the axial symmetric case), such that the two 
nucleons are subject to different densities. This should be 
contrasted with the case where the in-medium NN inter- 
action is considered in the azimuthal direction, namely 
where both nucleons are at the same distance from the 
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center of the nucleus and therefore subject to the same 
density. This is not a serious limitation since the latter 
interaction ought to be subdominant to the interaction in 
the radial direction, as the medium modifications due to 
the changes in the profile functions are known to be small, 
whereas the changes resulting from the density depen- 
dence of the mass functionals are more pronounced j^] . 



II. THE MODEL 

Unfortunately, the NN potential of the original Skyrme 
model does not have a central attraction at intermedi- 
ate distances. However it has been shown that this long 
standing problem can be solved e.g. by the incorpora- 
tion of the iterated two-pion exchan ge pol Isif into the 
Skyrme model or by a reformulation [S^ of the original 
Skyrme lagrangian on the basis of scale invariance and 
the conformal anomaly ■ Such a scale-invariant la- 

grangian was also used to describe nucleons in the nuclear 
matter |35ll36l | in the mean-field approximation (MFA). 



as X — Rk + where r is the distance vector relative to 
the origin of the skyrmion.^ We have chosen this model 
as a representative for a larger class of Skyrme-type mod- 
els that allow for a qualitative description of the central 
attraction between two nucleons. It is not our aim to 
construct a fine-tuned version that produces quantitative 
fits of the single nucleon properties or the two-nucleon 
potentials. Rather, by keeping the involved terms to the 
bare minimum for this case (i.e., the non- linear sigma- 
model kinetic term, the stabilizing fourth-order deriva- 
tive term, the symmetry-breaking pion mass term, and 
the intermediate-attraction generating dilaton term and 
couplings) this model should be general enough to pre- 
dict qualitatively those results which are generic and also 
hold for more complicated Skyrme-type lagrangians in- 
volving e.g. vector-mesons and other hadronic fields. 

The dependence on the nuclear density p{x) is included 
in the in- medium coefficients as{x) and ap{x), 



_ -, _ 47rcop(x)/77 
"^^"^^ l+5(,4^cop(f)A7' 
as{x) = 1 - A'n-ribopix)/ml . 



(2) 



A. Lagrangian and baryon number one soliton 

We start here with the model of Refs. js^l that proposes 
a scale-invariant version of the in-medium Skyrme model 
given, in the static case, by the following lagrangian: 
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-%^x'a.(x)Tr (2 -[/-[/+) . (1) 

Here U = exp{2ir-7f/i^^} parameterizes the pseudoscalar 
isotriplet of pion fields tt, t are the usual Pauli matrixes, 
X = e"^"^/^" is given in terms of the scalar-isoscalar dila- 
ton field (T, -Ftt is the weak pion decay constant, e is the 
parameter of the stabilizing Skyrme term, is the pion 
mass, £ = 16/29, and C* is the gluon-condensate pa- 
rameter.^ We associate the dilaton with quarkonium. 
Its mass rria- « 600 MeV is compatible with the well 
known attraction in the central NN potential [s^. The 
skyrmion is assumed to be located at a position Rk from 
the center of the nucleus, such that the total spatial vec- 
tor measured relative to the center of the nucleus is given 



Here rj = 1 + m^/mTv ~ 1-14 is a kinematical factor, 
and rriTv = 938 MeV is the mass of the nucleon. More- 
over, bo = — 0.024 and cq = 0.21 m""^ are empirical 
parameters that can be taken from the analyses of pio- 
nic atoms and low-energy pion-nucleus scattering data, 
and (7q = 1/3 is the Lorentz-Lorenz factor that takes into 
account the short-range correlations (38i] . 

Numerical calculations show that the ground state 
B=l solution of the Euler-Lagrange equations of (QJ in 
free space {p = 0, as = ap = 1) is spherically sym- 
metric. However, when the soliton is embedded into a 
finite nucleus, it may alter its shape, since the spherical 
symmetric configuration may not correspond to a mini- 
mum in energy or mass any longer. In such a way it is 
possible to study the modification of the shape of the in- 
medium skyrmion, namely by a minimization procedure 
of its mass functional. 

In deriving the mass functional of the single skyrmion 
we follow the scheme presented in Ref. by considering 
an axially symmetric configuration for a single skyrmion 
located at some distance from the center of a finite spher- 
ical nucleus. This implies that the chiral profile function, 
which parameterizes the modulus of the pion, is axially 
symmetric, i.e. F(f) = F{r,6), while the pion direction 
is governed by the polar-angle profile function Q{9): 

m = ^F(r,0)7V(e(0),^), 

N = {sin 6(61) cos (p, sin 9(6') sin ^, cos 6(6*)}. (3) 



The asterisk indicates that the gluon condensate should be con- 
sidered as renormalized one in the nuclear medium. Note that 
in general e = 8Nf/{llNc — "^Nf). Furthermore, we assume 

Fit = Fa-. 



^ This corresponds to U = U(x — Rk) = U{r) and x = x{^~^k) = 
x{?0 ^'i- Q- See Rcf. 5] for further details on the geometry 
of a skyrmion (without dilaton) inside a finite nucleus. 
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We further assume that the dilaton field also has an ax- 
ially symmetric configuration, i.e. cr(r) — a(r, 9) and 
therefore x(r) = xir^^)- For example, the mass func- 
tional for single skyrmion {B=l) located at a distance 
Rk =\Rk\ from the center of the nucleus and with a 
spatial vector fk measured relative to the origin of the 
skyrmion,^ is given as 



M{Rk) = 2^ y dr^rl j dOkSinOk 




l?2 



X |-^0(F,e;rfc,0fc) x'ap(x(rfc,0fc;i?fc)) 
1 



p2 2 

+ ^^(l-cosF) x'«s {xirk^eu^Rk)) 



C. 



(4) 



dilaton field: 



X = exp <^ - — (1 + 771 cos6'fc 



772 COS^ dk 



+ ...) 



(6) 

Here rg, 7^, Xd, fa and rji are variational parameters. 
As a result of the minimization procedure one gets a set of 
variational parameters corresponding to specified values 
of the distance and the other input parameters defined 
in subsection III CI and section ITTll 

In this context we should reiterate that we did not 
choose the model such that it reproduces in a fine-tuned 
way quantitatively the single nucleon properties or the 
two-nucleon interaction in free space, but that it rather 
serves as generic representative of Skyrme-type models by 
describing qualitatively the modifications when the single 
skyrmion or the two-skyrmion system is embedded in a 
finite nucleus. 



B. Two nucleons in a baryon-rich environment and 
the residual NN interaction 
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Here the subscripts of the functions F^^. , Fg^ , Qg^ , Xrk -i 
and X6k denote the corresponding partial derivatives (e.g. 
XBk = dx{rk,0k)/d9k of the function x = x{rk,Ok))- 

In Ref. the minimization of the mass functional Q 
was performed without dilatons, i.e. a{fk)=0, and the 
pertinent in-medium skyrmion properties were discussed. 
In particular, it was shown that the medium effects cause 
a deformation of the nucleon inside finite nuclei. We fol- 
low Ref. Q in the minimization procedure for the func- 
tional But in addition to the parameterizations of 
the chiral and polar-angle profile functions (5| 



F = 2 arctan 1 J [1 + 71 cos 6*^ + 72 cos^ 
e = Ok + Si sin 20k + S2 sin 46 k + S3 sin 66' a; - 



(5) 



we also introduce the following parameterization of the 



Note that here r 



rj, where k labels the skyrmion. Thus x 



= rk+Rk and \x\ = y/ Rl + rl + 2Rkrk cos Ok 
where and 9k are the modulus and polar angle of fk 



We now turn to the system of two interacting in- 
medium nucleons, e.g. to the modification of NN poten- 
tial Vnn inside nuclei. As mentioned in the introduc- 
tion, we only consider the modification of the NN poten- 
tial in the radial direction, since the modification in the 
azimuthal direction is expected to be subleading. The 
geometry of this axially symmetric case is presented in 

Figni 




FIG. 1: The axially symmetric two-skyrmion system in a fi- 
nite spherical nucleus. Here O is the center of the nucleus, 
O' is the geometrical center of the two-skyrmion system, Oi 
and O2 are the centers of the first and second skyrmions, re- 
spectively. D is the distance between the center of the nucleus 
and the one of the two-skyrmion system, and R is the distance 
between skyrmions. 

Unfortunately, in any version of the Skyrme la- 
grangian, it is very difficult to derive the binding energies 
as well as the form factors of even light nuclei directly 
from the nucleon- nucleon potential Vat at. For simplicity, 
in setting up our formalism, we assume that 
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• each nucleon is bound inside the finite nucleus by 
a phcnonicnological averaged potential normally 
used in shell models in MFA; 

• the nucleon distribution, i.e. the nuclear density, is 
given by a phcnomcnological formula; 

• the primary two-nucleon potential in free space 
is described by the lagrangian of the Skyrme 
model ^ with ap — as — 1- 

Although for practical calculations of nuclear properties 
in MFA the explicit expression of Vnn is not needed, this 
does not mean that the nucleons do not interact with each 
other. In fact, it is believed that the sum of all two-body 
- as well as possible three-body - interactions can be re- 
expressed as an effective averaged potential, centered at 
the origin of the nucleus. Due to the MFA ideology, every 
nucleon behaves as if it "feels" only this effective poten- 
tial (spherical- well potential, oscillator, etc.). Yet resid- 
ual interactions between nucleons may still exist. These 
interactions ought to be important for the description of 
deformations and stabilities of the nuclei. 

In the following we will study the residual part of NN 
interactions (in radial direction) embedded in the sur- 
rounding nuclear environment. We assume that the pri- 
mary two-nucleon potential in free space is given by the 
Skyrme model, including dilatons, as described by the la- 
grangian 1^ with density p — 0. We then "plug" a pair 
of nucleons into a finite nucleus, or more exactly, we con- 
sider them as part of the nucleus. The point is that each 
of the nucleons can be deformed for the following two 
reasons, which affect the in-medium potential as well: 

1) the presence of the other nucleon or, in other words, 
the presence of the NN interaction; 

2) the presence of the remaining A— 2 nucleons or, in 
other words, the presence of medium effects in the 
nucleus. 

The first case was investigated in Refs. 0,113, where the 
modifications of the nucleons and Vnn properties in free 
space were recorded under possible deformations of the 
nucleons during their mutual approach. For that purpose 
the authors introduced a couple of deformation param- 
eters and obtained a modified central potential V^jv 
minimizing it at each internucleon distance R. The sec- 
ond case, namely the effects of the nuclear medium, is 
the subject of this article. In pri ncip al, we could have 
followed the procedure of Refs. |1,|33, which would be: 

1) to introduce deformation parameters in la- 
grangian ^ using an axially symmetric ansatz; 

2) to construct V^^v from the lagrangian in the 
product approximation; 

3) to determine the deformation parameters by mini- 
mizing V^j^ at each R. 



However such a procedure would be unnecessarily com- 
plicated, especially for the case of lagrangian with 
dilatons. Moreover, it that case one would have to fight 
complications resulting from double counting the defor- 
mation effects caused by NN interactions and by medium 
modifications, respectively. Instead, for simplicity and 
also for isolating the medium modifications from the 
vacuum NN interaction effects, we minimize the mass 
functional 10} for each of the two in-medium skyrmions 
with respect to the corresponding deformation parame- 
ters, and then we insert these deformed skyrmion profiles 
into the product ansatz in order to generate the NN cen- 
tral potential. 

In this way, by inserting two skyrmions in the nuclear 
background of the A— 2 remaining nucleons, one can de- 
fine a static skyrmion-skyrmion residual potential 

Vssil\Ri + R2\, \R2-Ri\) 

= M[Uix_^)U{x_^),x{ri)x{r2);PA-2i\m 

=ri =f2 

-2MiT^-5V^I^{D,R). (7) 

Here M[nJ^^ U{x — Ri), x{x — Ri); pA-2i\x\)] is the static 
mass functional of a skyrmion-product ansatz of baryon 
number B — 2 with the skyrmions centered (axial- 
symmetrically) at Ri = OOi (with O and Oi as de- 
fined in FigQ in the renormalized background density 
Pyi-2(|5'|)- Thus the potential Q is defined in analogy 
with the well known product-ansatz prescription of the 
dilaton-free case 

Ub=2[x~ D,ltj ^U{x~ Ri)U{x~ R2), (8) 

where D = {Ri + R2)/2 and R = R2 — Ri- The remain- 
ing quantity 6Vgl^ denotes the contribution of the B=2 
system to the mean-field potential. It is defined as 

2 

i=l 

(9) 

Note that the corresponding quantity 6 Ms = —SV^i^ 
can be interpreted as the part of the energy of skyrmion 
that is transformed to the mean-field energy of pions 
inside the nucleus and to the interaction energy of the 
skyrmion with the background field. It is this back- 
ground pionic field that generates the attractive mean 
field potential binding nucleons to the nucleus. Clearly 
the potential O becomes the free NN potential when 
the background density of the remaining nucleons goes 
to zero"*. 



* If the background field vanishes, pA-2 — > 0, the mean-field con- 
tribution also vanishes, SVgiJ!^ ~* 0- 
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Finally, inserting © into Q, one gets the following 
form of the skyrmion-skyrmion residual potential 

Vss{D,R) = M[U{n)U{f2),x{ri)x{r2);PA-2{\m 
-M[U{ri),x{ri);PA-2{\m 
~M[U{r2),x{r2hPA-2{\m- (10) 

From the comparison of Eqs. (0), ©, and it should 
become clear that the potential (|10|) is simply the differ- 
ence between the separation energy of a pair of skyrmions 
{B=2) inside the nucleus and the sum of the separation 
energies of each individual skyrmion (B—l) in the same 
nucleus. This is well known in nuclear physics as a major 
contribution to the residual NN potential in the nuclear 
medium. Our aim is the evaluation of the spin-isospin 
independent part of this potential, potential. Note that 
in this case the skyrmion-skyrmion interaction and the 
nucleon-nucleon interaction coincide. 

As mentioned before, we only consider the geometry 
when the center of the (spherical) nucleus and centers 
of the two skyrmions fall onto the same axis, which is 
denoted by z as illustrated in Fig^ Then the V^n po- 
tential depends only on two variables: 1) the distance be- 
tween the centers of the two skyrmions, R; 2) the distance 
between the geometrical center of the two-skyrmion sys- 
tem and the center of the nucleus, D. The spatial vector 
r is measured relative to the two-skyrmion center, such 
that here x — D + f. 

To isolate the spin-isospin independent part of the 
NN potential one can use the standard projection pro- 
cedure [40j |. which corresponds to the following form in 
free space 



Vnn{R, C) 



NN 



-f -Tr(Co-,;C+aj) pi?,;i?j-(5,y) 
xV^AR)+5.,V^-^^{R)\, (11) 



where C = uq -\- if ■ u (with + — 1) \s the relative 
orientation matrix of the nucleons in internal space, Ci j- 
are the Pauli matrixes, and R is the unit vector along the 
line joining the centers of the two nucleons. Here V^j^, 
V^-^ and V^]v are the central, tensor and spin-spin parts 
of the NN potential, respectively. 

Finally, taking into account the projection proce- 
dure 40] and using the expression for skyrmion-skyrmion 
potential H1U|) , we predict the following expression for the 
spin-isospin-independent part of the residual NN poten- 
tial (in radial direction) inside a finite spherical nucleus: 

OO TT 

y^'^(D,i?) = 2^ Jdrr'Jdesme{V,^ + V,' + V^sB + V:} 



(12) 



with 



^2 = ^-faj,,A-2{x{r,9;D))[xlxU(t>i+h) 

- (X?01 +X2'^2)], 



-^{M2-^) 



v: 



XSB 



-2 {x[r 



\0-,D)) [x?X^(l-cosFiCOsF2) 
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-X?(l-cosFi)-X2(l-cosF2)] , 
iAxi ■ X2) - V'(Xi) - V'(X2) 



'a,A- 2 
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-^XiX2 



sin( 



/ 

) ^ N / X2,e2 Xifii \ 

'l-6'2 Xl,ri X2,r2 

\ r2 n J 



$ = 



+ COs(0i-02) (xi,r.,X2,r. + ^^Ml^ 

h-02) (F,.r,^~^F2,r^ 

\ r2 n J 



sin 



cos( 



^1 — ^2) f -Fl.ri^2.r2 H h2j—:£l\ 

\ "^1^2 ) _ 




- sin( 



01 - 92)Fi^ri I 62,9 



-^Fi sin^ F2 ( sir? 61 sin^ 62 



sin 



cos2(0i-02)e?,e^e2_, 




Here Fi = Fi{n,e^), Xi = Xii^i^Oi), 9, = e,(6'i) with 
i — 1,2 are the chiral, dilaton, and polar- angle profile 
functions of skyrmion 1 or 2, respectively, with — 
ri{r,e;R) = yf {K^ / A) + ± Rr cosO, 9, = 9i{r,9;R) = 
arcsin(sin0r/ri(r, 0; i?)), r, 9 and x — x{r,9;D) = 
y/D^ +r^ + 2Dr cos 9 as defined in Fig.^ Moreover, 
we abbreviate the partial derivative of the above defined 
functions G {Fi,Xi,Qi} as fi^y = dfi/dy. The func- 
tions (j)i are defined as <f>i = 4>{Fi,Qi;ri,9i) with (j) (and 
■0) as in Eq. (@J. 



C. Density dependence and input parameters 

We shall use a phenomenological parameterization of 
the density of a spherical nucleus as presented in Ref. ^3 
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and therefore write the renormahzed density of the A — B 
nuclear background environment as 



Pa-b{x) = 



Pa-b{x) 



A~B 
A 

X exp 
A-B 



7r3/2,,3 



1 



A-2 



iiA< 20, 



Pq 



A J l + exp{ (a; -i?')/a} 



if A > 40. 



where ro = 
0.58 fm, R' 



1.635 fm as for ^^C. Furthermore, a = 
: 1.2A^/^fm, and po ~ 0.5m^ is the nor- 
mal nuclear matter density. 

The input parameters of the skyrmion sector are cho- 
sen as =186 MeV and e = 27r. Although this standard 
set of input parameters gives an overestimated value for 
the mass of nucleon even in free space, we shall not opti- 
mize it, since we are only interested in the influence of the 
medium on the NN potential. The sole input parameter 
in the dilaton sector is the gluon condensate parameter 
C*. Its in- medium renormalization has not been clari- 
fied in the literature. In the present model this quantity 
can be expressed in terms of the in-medium sigma-meson 
mass TO* as IsIIstII: 



C. 



A-B 



g.A-1 



2(4 -e) 



(13) 



Unfortunately, there is also no information about the 
renormalization of to ^ w ithin the present approach. Var- 
ious approaches [s^, jl^l show that m* has a linear 
density dependence. In view of this, we shall use the 
parameterization 



1_ 0.12^^^') TO. (14) 
Po 



that was obtained in Ref. 43] within the quark-meson 
coupling model. ^ The free-space mass of the sigma meson 
is taken as to. = 550 MeV, which corresponds to a free- 
space value of Cg = (260 MeV)''. 



III. RESULTS AND DISCUSSION 

The results of the minimization procedure for a single 
skyrmion (i? = 1) in the nucleus with an yl — 2 back- 
ground density is presented in Tabled We have not 
listed all variational parameters, as some of them are 
rather small, namely 7^, Si, rji for i > 3. These param- 
eters, which describe higher multipole deviations of the 
skyrmion shape from the spherical form, can be neglected 



^ Since in finite nuclei the density is coordinate-dependent, 
will also acquire a coordinate dependence. 



Rk[im] 





±0.649 ±1.298 


±1.947 


±2.595 


rs [fm] 


0.512 


0.509 0.495 


0.467 


0.435 


71 


±0.020 tO.034 =f0.105 t0.107 


72 





-0.028 -0.047 


0.001 


0.037 


105i 





-0.056 -0.097 


-0.018 


0.047 


WS2 


O 


-0.002 


-0.003 


-0.001 


Xd 


0.634 


0.629 0.592 


0.655 


0.766 


rd[fm] 


0.843 


0.841 0.835 


0.808 


0.776 


Vi 


±0.055 ±0.008 tO.060 t0.054 


V2 





-0.015 -0.032 





0.017 




Rk[im] 





±1.705 ±3.409 


±5.114 


±6.819 


rs [fm] 


0.597 


0.576 0.477 


0.407 


0.398 


71 


=f0.053 tO.145 t0.035 =f0.002 


72 





-0.011 0.030 


0.017 


0.001 


105i 





-0.026 0.021 


0.027 


0.002 


WS2 





-0.001 








Xd 


0.397 


0.423 0.650 


0.854 


0.876 


rd[fm] 


0.983 


0.951 0.818 


0.749 


0.741 


Vi 


=f0.020 tO.076 tO.010 tO.OOI 


V2 





-0.005 0.014 


0.004 






TABLE I: Deformation parameters at various separations 
Rk of the center of a single skyrmion from the center of the 
nucleus. The parameters rs, 7i, Si are for the chiral sector, 
while Xd, fd, r]i are for the dilaton sector. They are obtained 
by a direct minimization of the mass functional in |0J. Note 
that all coefficients 7^, Si, rji for i > 3 are negligible and 
are not presented here. Non-zero free-space values of some 
parameters are: rs = 0.398 fm, Xd = 0.877, rd = 0.740 fm. 



without loss of accuracy. A discussion about the param- 
eters of the chiral sector can be found in Ref. |^] , whereas 
the ones of the dilaton sector are new. Both sets are con- 
sistent with the fact that the nucleon inside the nucleus 
acquires an intrinsic quadrupole moment unless the nu- 
cleon is located directly at the center. Furthermore note 
that the increase of the radial parameter rs of the chi- 
ral sector and of the radial parameter of the dilaton 
sector with increasing density is in qualitative agreement 
with the overall picture of the swelling of the nucleon 
inside the nucleus. The behavior of the depth of the 
dilaton field parameterized by Xd is more involved, but 
at least for the case of the heavier nucleus the dilaton 
region of the in-medium nucleon becomes shallower with 
increasing density. At first sight this is counter-intuitive. 
However, note that because of the swelling of the nucleon 
the total region where the "effective pion decay parame- 
ter" xF-rr is reduced relatively to its vacuum value Ftt has 
increased nevertheless. 

The results of the calculations for the residual NN in- 
teractions are summarized in Figs. [21 and |21 The spin- 
isospin-independent central part V^j^ of the residual NN 
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FIG. 2: Spin-isospin-independent central part (in radial di- 
rection) of the residual NN potential in ^'^C. The solid line 
corresponds to the potential in the free case, the dashed line 
corresponds to Z) = 3 fm, the dotted line corresponds to 
D = 2 fm, and the stars correspond to D = fm. 
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FIG. 3: Spin-isospin-independent central part (in radial di- 
rection) of the residual NN potential in *''Ca. The solid line 
corresponds to the potential in the free case, the dashed line 
corresponds to Z) = 4 fm, the dotted line corresponds to 
D = 3 fm, and the stars correspond to D = fm. 



potential in ^^C is shown in Fig.|21 Let us emphasize 
again that our model has not been fine-tuned to repro- 
duce the spin-isospin independent NN potential in free 
space quantitatively. But the qualitative changes are ev- 
ident: One can see that short-range repulsive contribu- 
tion of the potential is increased and the intermediate 
attraction is decreased in the nuclear medium. The same 



behaviour of is observed for skyrmions embedded 

in ^°Ca, as shown in Fig.|21 

It is known that the residual NN potential becomes re- 
pulsive when neighboring nucleons overlap. It is further- 
more known that the sizes of the nucleons increase in the 
medium, see also the parameters rs and of Table 
Therefore the overlap-regions increase with density. This 
leads to a build-up of the repulsive part of the poten- 
tial at high densities, i.e. when the two-nucleon system 
is near the center of the nucleus. 

From previous studies Q we know that the nucleoli 
mass is the smallest near the center of nucleus and in- 
creases when the nucleon is moved to the surface. Due 
to this phenomenon, nucleons should actually collapse to 
the center of the nucleus. In other words, due to the 
attractive mean-field potential, nucleons should move to 
the center. On the other hand, the nucleon concentration 
near the center of the nucleus does, as shown in Figs.[21 
and 131 increase the repulsive residual potentials between 
nucleons such that nucleons are expelled from the center. 
Consequently, an equilibrium state arises and therefore 
saturation of the nuclear matter density results. In this 
case nucleons stop their radial motion towards the cen- 
ter of the nucleus, but their motion in a shell with given 
radius continues. A further inclusion of the Pauli mech- 
anism in this picture could finally provide for the shell 
description of finite nuclei. 

Su et al. obtained similar results to ours in the frame- 
work of a chiral a-to model under the imaginary-time 
Green's function method • They found that the po- 
tential well of the nucleon-nucleon interaction becomes 
shallow as the temperature increases. The same be- 
haviour of the NN potential was seen in quantum hadro- 
dynamics studies and in the framework of TFD ^2§\ . 
The latter established, in addition, the alternate roles of 
density and temperature effects. 



IV. SUMMARY AND OUTLOOK 

We have considered the central part of the residual 
NN potential inside finite nuclei in the radial direction. 
In the interior region of the heavy nucleus the residual 
NN potential is strongly repulsive and compensates the 
mean-field attractive potential in a such way that the 
nuclear density saturates. These results are consistent 
with other studies applying different methods 0, 
lisf . However, one should remember that the Skyrme- 
type models are the only models that are formulated in 
terms of hadronic (non-quarkish) degrees of freedom and 
that can simultaneously describe both single-baryon as 
well as multi-baryon properties. In the framework of such 
a class of models our results are indeed new. 

It would be interesting to investigate the role of the 
residual interactions in the formation of nuclear mat- 
ter itself, i.e. to study their role in the formation of 
the mean-field potential. This requires the additional 
study of the in-medium NN potential in connection with 
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nuclear matter properties. For example, for some rela- 
tive orientations of nucleons in the internal space (see 
Eq. (fTT|l ) the tensor part of the potential V^f^ makes a 
negative contribution to the total NN potential. Such an 
attraction turned out to be sufficient for the construction 
of crystalline |^ or condensate (46|] states of nuclear mat- 
ter. Note, however, that medium modifications of the 
NN potential were neglected. Previous calculations 
showed that the tensor part of the residual potential does 
decrease in the nuclear medium. Thus the Skyrme model 
allows to combine the tensor and the central potentials 
into a full potential which is shallow due to the medium 
influence. The question is whether the changes of the 
in-medium NN potential are sufficient to induce a break- 
down of the crystalline structure which seems not to exist 
in nature. Of course the medium modifications of the NN 
potential might only be one of many agents that could 



induce this breakdown. For instance, quantum fluctua- 
tions p6L and Fermi- motion effects should be of im- 
portance |48j| . It is well known that the nuclear matter 
binding results from a strong cancellation of an attractive 
(binding) potential term and a repulsive kinetic term. It 
remains to be seen whether the condensed state of nu- 
clear matter at ordinary densities can be achieved within 
the in-medium modified Skyrme model. 

This research is part of the EU Integrated Infra- 
structure Initiative Hadron Physics Project under con- 
tract number RII3-CT-2004-506078 and is supported in 
part by DFG (SFB/TR 16) and by the Forschungszen- 
trum Jiilich under contract No. 41445400 (COSY-067). 
U.T.Y. would like to thank A.K. Nasirov for useful corre- 
spondence and acknowledges support by Pusan National 
University (Post-Doc 2004 program). 



[1] A. Jackson, A.D. Jackson and V. Pasquier, Nucl. Phys. 

A 432, 567 (1985). 
[2] U.B. Kaulfufi and U.-G. Meifiner, Phys. Rev. D 31, 3024 

(1985) . 

[3] C. Hajduk, B. Scliwesinger, Nucl. Phys. A 453, 620 

(1986) . 

[4] T. Otofuji, S. Saito, M. Yasuno, H. Kanada and T. Kuri- 

hara. Prog. Theor. Phys. 78 527 (1987). 
[5] U.T. Yakhshiev, M.M. Musakhanov, A.M. Rakhimov, 

U.-G. Meifiner, A. Wirzba, Nucl. Phys. A 700, 403 

(2002). 

[6] U.T. Yakhshiev, U.-G. Meifiner, A. Wirzba, Eur. Phys. 

J. A 16, 569 (2003). 
[7] A.M. Rakhimov, M.M. Musakhanov, F.C. Khanna and 

U.T. Yakhshiev, Phys. Rev. C 58, 1738 (1998). 
[8] I.R. Klebanov, Nucl. Phys. B 262, 133 (1985). 
[9] N.K. Glendenning, Phys. Rev. C 34, 1072 (1986). 
[10] E. Wiist, G.E. Brown, A.D. Jackson, Nucl. Phys. A 468, 

450 (1987). 

[11] A.S. Goldhaber and N.S. Manton, Phys. Lett. B 198, 231 

(1987) . 

[12] M. Kugler and S. Shtrikman, Phys. Lett. B 208, 491 

(1988) ; M. Kugler, S. Shtrikman, Phys. Rev. D 40, 3421 

(1989) . 

[13] A.D. Jackson and J.J.M. Verbaarschot, Nucl. Phys. A 
484, 419 (1988). 

[14] L. Castillejo, P.S.J. Jones, A.D. Jackson, J.J.M. Ver- 
baarschot, A. Jackson, Nucl. Phys. A 501, 801 (1989). 

[15] H. Forkel, A.D. Jackson, M. Rho, C. Weiss, A. Wirzba 
and H. Bang, Nucl. Phys. A 504, 818 (1989). 

[16] W.K. Baskerville, Nucl. Phys. A 596, 611 (1996). 

[17] N.S. Manton, P.J. Ruback, Phys. Lett. B 181, 137 
(1986). 

[18] N.S. Manton, Commun. Math. Phys. Ill, 469 (1987). 

[19] A.D. Jackson, A. Wirzba, L. Castillejo, Phys. Lett. B 
198, 315 (1987); A.D. Jackson, A. Wirzba, L. Castillejo, 
Nucl. Phys. A 486, 634 (1988); A.D. Jackson, C. Weiss, 
A. Wirzba and A. Lande, Nucl. Phys. A 494, 523 (1989). 

[20] A.D. Jackson, N.S. Manton and A. Wirzba, Nucl. Phys. 
A 495, 499 (1989). 

[21] A. Wirzba and H. Bang, Nucl. Phys. A 515, 571 (1990). 



[22] A.D. Jackson, C. Weiss, A. Wirzba, Nucl. Phys. A 529, 
741 (1991). 

[23] A. Wirzba, in: G. Holzwarth (Ed.), Baryons as 
Skyrme Solitons, World Scientific, Singapore, 1993, 
arXiv:hep-ph/9211295 

[24] B.D. Serot, J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 

[25] S. Gao, Y.-J. Zhang and R.-K. Su, Nucl. Phys. A 593, 
362 (1995). 

[26] Y.-J. Zhang, S. Gao and R.-K. Su, Phys. Rev. C 56, 3336 
(1997). 

[27] L. Alvarez-Ruso, P. Fernandez de Cordoba and E. Oset, 

Nucl. Phys. A 606, 407 (1996). 
[28] A.M. Rakhimov, U.T. Yakhshiev and F.C. Khanna, 

Phys. Rev. C 61, 024907 (2000). 
[29] U.T. Yakhshiev, A.W. Thomas and F.C. Khanna, Phys. 

Rev. C 68, 048201 (2003). 
[30] N. Kaiser, U.-G. Meifiner, Phys. Lett. B 233, 457 (1989). 
[31] N. Kaiser, U.-G. Meifiner, Nucl. Phys. A 506, 417 (1990). 
[32] M.M. Musakhanov, A.M. Rakhimov, Mod. Phys. Lett. A 

10, 2297 (1995). 
[33] H. Gomm, P. Jain, R. Johnson and J. Schechter, Phys. 

Rev. D 33, 3476 (1986); A.A. Migdal and M.A. Shifman, 

Phys. Lett. B 114, 445 (1982). 
[34] V.A. Andrianov, V.Yu. Novozhilov, Phys. Lett. B 202, 

580 (1988). 

[35] G.E. Brown, M. Rho, Phys. Rev. Lett. 66, 2720 (1991); 

M. Rho, Phys. Rep. 240, 1 (1994). 
[36] G.E. Brown, M. Rho, Phys. Rep. 269, 333 (1996). 
[37] M.M. Musakhanov, A.M. Rakhimov, U.T. Yakhshiev, 

Z. Kanokov, Phys. Atom. Nucl. 62, 1845 (1999); 

A.M. Rakhimov, F.C. Khanna, U.T. Yakhshiev, 

M.M. Musakhanov, Nucl. Phys. A 643, 383 (1998). 
[38] T.E. Ericson and W. Weise, Pions and Nuclei (Claredon, 

Oxford, 1988). 

[39] A.M. Rakhimov, T. Okazaki, M.M. Musakhanov and 
F.C. Khanna, Phys. Lett. B 378, 12 (1996). 

[40] E.M. Nyman and D.-O. Riska, Phys. Scr. 34, 533 (1986). 

[41] A.I. Akhiezer, A.G. Sitenko, V.K. Tartakovski, Nuclear 
electrodynamics (Springer Verlag, 1994). 

[42] U.-G. Meifiner, Nucl. Phys. A 503, 801 (1989). 

[43] K. Saito, K. Tsushima and A.W. Thomas, Phys. Rev. C 



9 



55, 2637 (1997). 
[44] R.-K. Su, S.-J. Yang, S. Gao, P.K.N. Yu, J. Phys. G 20, 
1757 (1994). 

[45] S. Gao, Y.-J. Zhang, R.-K. Su, Phys. Rev. C 52, 380 
(1995). 

[46] D. Diakonov and A.D. Mirhn, Sov. J. Nucl. Phys. 47, 21 



(1988) [Yad. Fiz. 47, 662 (1988)]. 
[47] T.S. Walhout, Nucl. Phys. A 484, 397 (1988); T.S. Wal- 

hout, Nucl. Phys. A 519, 816 (1990). 
[48] T.D. Cohen, Nucl. Phys. A 495, 545 (1989). 



